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Abstract 

We study the two-body and three-body charmful baryonic B decays: B° —* A c p and B — > 
A c pir(p). The factorizable VF-exchange contribution to 5° — ► A c p is negligible. Applying the 
bag model to evaluate the baryon-to-baryon weak transition matrix element, we find B(B° — ► 
A-cp) ^ 1.1 x 10 _5 |5 r B0 pS + /6| 2 with g B o p Y;+ being a strong coupling for the decay — ► B°p and 
hence the predicted branching ratio is well below the current experimental limit. The factorizable 
contributions to B~ — > A c pit~ can account for the observed branching ratio of order 6 x 10 -4 . The 
branching ratio of B~ — * A c pp~ is larger than that of B~ — * A c pir~ by a factor of about 2.6 . We 
explain why the three-body charmful baryonic B decay has a larger rate than the two-body one, 
contrary to the case of mesonic B decays. 
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I. INTRODUCTION 



Inspired by the claim of the observation of the decay modes ppir^ and pp-n^-n^ in B 
decays by ARGUS ||T| in the late 1980s, baryonic B decays were studied extensively around 
the early 1990s with the focus on the two-body decay modes, e.g. B — > pp, AA. Up 

to now, none of the two-body baryonic B decays have been observed. Indeed, most of the 
earlier predictions based on the pole model or QCD sum rule or the diquark model are too 
large compared to experiment fll^JTol (see Table I). 



TABLE I. Predictions of branching ratios for some two-body baryonic B decays in various 
models. We have normalized the branching ratios to [l^/V^I = 0.085. The predictions given in 
[11] are carried out in two different quark-pair-creation models: local and nonlocal. Experimental 
limits are taken from [14,15]. 
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[11] 

non-local local 


experiment 




■* A c p 






4 x 10" 4 


8.5 x 10~ 4 


1.1 x 10" 3 


1.7 x 10~ 3 1.9 x 10~ 3 


< 2.1 x 10" 4 


B°- 


^pp 




4.2 x 10" 6 


1.2 x 10~ 6 


5.9 x 10~ 6 


7.0 x 10~ 6 


2.9 x 10~ 6 2.7 x 10~ 5 


< 1.6 x 10" 6 


B°- 


-> AA 








3.0 x 10~ 6 




1.2 x 10~ 6 2.3 x 10~ 5 


< 2.3 x 10~ 6 


B- - 


-> A— 


~P 


1.5 x 10~ 4 


2.9 x 10~ 7 




3.2 x 10~ 4 


7.2 x 10~ 7 8.7 x 10~ 6 


< 1.5 x 10~ 4 



In order to understand why the momentum spectrum of produced A c in inclusive B 
decays is soft and why the two-body decay modes, e.g. B — > {A c , S c }{p, A}, have not been 
observed, Dunietz |T(| argued that a straightforward Dalitz plot for the dominant b — > cud 
transition predicts the cd invariant mass to be very large. The very massive cdq objects 
would be usually seen as A c nn(n > 1) if the cd forms a charmed baryon. This explains the 
observed soft A c momentum spectrum and the non-observation of A c p decay. Since the very 
massive cdq could also be seen as D^NX, the baryonic processes B — ► D^NN'X would 
be likely sizable. Indeed, CLEO has recently reported the observation of B° — > D*~pn at the 
10~ 3 level and B° — > D*~ppix + at the 10~ 4 level [|17|. Theoretically, the three-body decay 
modes B — > D*~NN and B° — > p^(ji~)pn have been recently studied in |18|Jl9| . 

A similar observation has been made by Hou and Soni f2~0fl . They pointed out that the 
smallness of the two-body baryonic decay B £>i£> 2 has to do with the large energy release. 
They conjectured that in order to have larger baryonic B decays, one has to reduce the 
energy release and at the same time allow for baryonic ingredients to be present in the final 
state. Under this argument, the three-body decay, for example B — > ppn, will dominate 
over the two-body mode B — »• pp since the ejected p meson in the former decay carries 
away much energy and the configuration is more favorable for baryon production because 
of reduced energy release compared to the latter |19|. This is in contrast to the mesonic 
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B decays where two-body decay rates are usually comparable to the three-body modes. 
The large rate of B° — > D*~pn and B° — > D*~pp-K + observed by CLEO indicates that the 
decays B — ► baryons receive comparable contributions from B — > A c pX and B — > [DjiViV'X, 
where [D] denotes any charmed meson. By the same token, it is expected that for charmless 
baryonic B decays, B — > p(ir)BiB2 are the dominant modes induced by tree operators and 
B — > (71", rf, p)Bu s )B2, e.g. B — > pAp, are the leading modes induced by penguin diagrams. 
In this work we focus on charmful baryonic decays B —>■ A c pX. The experimental results 



are summarized as 21 



B[B- -> A c p7T-) = (6.2 ± 2.7) x 10" 4 , B{B° -> A c p7r ) < 5.9 x 10~ 4 , 
£(5° -> A c p7r + 7r") = (1.3 ±0.6) x 10" 3 , B(B~ ^ A c pir- ir°) < 3.12 x 10- 3 , (1.1) 

to a ^ . „ , 4 



together with the upper limit £>(£> — > A c p) < 2.1 x 10 _ . It is evident that the two-body 
mode is suppressed. Specifically, we shall study B° — > A c p and B~ — > A c p-K~ (p~) in detail 
in order to understand their underlying decay mechanism. It has been advocated that the 



B decay to A c p + 7r's is suppressed relative to A c p [10]. We shall see that this is not the case. 



The layout of the present paper is organized as follows. In Sec. II we first study the 
two-body charmful decay B — > A c p to update the prediction of its branching ratio. We 
then turn to the three-body decays B — > A c pn(p) in Sec. III. A detail of the MIT bag model 
for the evaluation of baryon-to-baryon weak transition matrix elements is presented in the 
Appendix. 

II. TWO-BODY CHARMFUL BARYONIC DECAY B° -> A c p 

We first study the two-body baryonic decay 5° — > A c p to update its prediction and 
understand why it is suppressed compared to three-body modes. To proceed, we first write 
down the relevant Hamiltonian 

H eS = %K6Kj[ci(p)0i(/x) + C2(jl)0 2 (jJL)] + h.c, (2.1) 

where 0\ = (cb)(du) and 2 = (cu)(db) with (g^) = <?i7^(l — 75)<?2- In order to ensure 
that the physical amplitude is renormalization scale and 75-scheme independent, we include 
vertex corrections to hadronic matrix elements. This amounts to modifying the Wilson 
coefficients by [[2211 : 



CM - of = cM + ^ U^ln^ + f T ) aQi), 



c 2 {p) c\ = c 2 (/x) 



+ S(^ ,r "V + fT )/-(")- < 2 - 2 » 



where the anomalous dimension matrix 7^ and the constant matrix f in the naive dimen- 



sional regularization and 't Hooft-Veltman schemes can be found in [22]. The superscript 
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T in Eq. (|2.2|) denotes a transpose of the matrix. Numerically we have c\ s = 1.168 and 
c| ff = —0.365 p2| . It should be stressed that cf^ and c| are renormalization scale and 



scheme independent For later purposes we write 



H eS = ^V cb V: d [c + + + c_0_] + h.c, (2.3) 

with 0± = Ox ± 2 and c± = |(cf ± cf ). 

The decay amplitude of 5° — > A c p consists of factorizable and nonfactorizable parts: 

A(B° -> A c p) = -> A c p) fact + -> A c p) nonfact , (2.4) 

with 

A(B^ - A c p) fact = ^V a V: d a 2 {A c p\(cu)\0){0\(db)\^), (2.5) 

where a 2 = c| ff + c\ s /N c . The short-distance factorizable contribution is nothing but the 

W-exchange diagram. This W^-exchange contribution has been estimated and is found to be 

very small and hence can be neglected [fjUl3|| . However, a direct evaluation of nonfactorizable 

contributions is very difficult. This is the case in particular for baryons, which being made 

out of three quarks, in contrast to two quarks for mesons, bring along several essential 

complications. In order to circumvent this difficulty, it is customary to assume that the 

nonfactorizable effect is dominated by the pole diagram with low-lying baryon intermediate 

states; that is, nonfactorizable s- and p-wave amplitudes are dominated by ~ low-lying 

baryon resonances and | + ground-state intermediate states, respectively PTJ| . For B° — > A c p, 

— o +(*) — 

we consider the strong-interaction process B — > E 6 y p followed by the weak transition 
S^*" 1 — > A c , where is a | baryon resonance (see Fig. 1). Considering the strong 



coupling 



i9B P Z b 1pV b l5lpp(l>B + gBpZ^^p^B, (2.6) 



the pole-diagram amplitude has the form 



where 



A(B° -> A c p) nonfact = u Ac {A + B^ 5 )v p , (2.7) 



*For the mesonic decay B — > MiM 2 with two mesons in the final state, two of the four quarks 
involving in the vertex diagrams will form an ejected meson. In this case, it is necessary to take 
into account the convolution with the ejected meson wave function. 
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FIG. 1. Some pole diagrams for 13° — > A c p where the symbol • denotes the weak vertex, (a) 
corresponds to nonfactorizable internal W emission, while (b) to the VF-exchange contribution. 

correspond to s-wave parity-violating (PV) and p-wave parity-conserving (PC) amplitudes, 
respectively, and 

(A c |K P ff C |^> = u Ac a EbAc u Eb , (A c \H p J\J:i + ) = m Ac 6 s . Ae « S *. (2.9) 

The main task is to evaluate the weak matrix elements and the strong coupling constants. 
We shall employ the MIT bag model J23| to evaluate the baryon matrix elements (see e.g. 



2~4j , |25|| for the method). Since the quark- model wave functions best resemble the hadronic 
states in the frame where both baryons are static, we thus adopt the static bag approximation 
for the calculation. Note that because the four-quark operator + is symmetric in color 
indices, it does not contribute to the baryon-baryon matrix element since the baryon-color 
wave function is totally antisymmetric. From Eq. ( |2.3| ) and the Appendix we obtain the PC 
matrix element 

a Sb A c = -^V cb V^^={X 1 + 3X 2 )(4tt), (2.10) 

where 

f R 

Xi= r 2 dr[u c (r)v u (r) - v c (r)u u (r))[u d (r)v b (r) - v d (r)u b {r)}, 
Jo 

f R 

X 2 = r 2 dr[u c (r)u u (r) + v c (r)v u (r)][u d (r)u b (r) + v d (r)v b (r)], (2.11) 
are four-quark overlap bag integrals (see the Appendix for notation). In principle, one can 



also follow |24| to tackle the low-lying negative-parity state in the bag model and evaluate 
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the PV matrix element &e*a c |] However, it is known that the bag model is less successful even 
for the physical non-charm and non-bottom | resonances [p3[| , not mentioning the charm 
or bottom | resonances. In short, we know very little about the ~ state. Therefore, we 
will not evaluate the PV matrix element &e*a c as its calculation in the bag model is much 
involved and is far more uncertain than the PC one |fJ3| . 

Using the bag wave functions given in the Appendix, we find, numerically, 



The decay rate of B - 
Y(B 



-1.49 x l(T 5 GeV 3 , 
— > B\B 2 is given by 

B X B 2 ) = 

+ \B 



X? = 1.81 x l(T 4 GeV 3 



(2.12) 




(m B + m 1 + m 2 ) 2 p 2 c 
(Ei + m 1 )(E 2 + m 2 )m 2 B 
2 [(E 1 + m 1 )(E 2 + m 2 )+p 2 ] 2 \ 
(Ei + mi)(E 2 + m 2 )m 2 B j 



(2.13) 



where p c is the cm. momentum, and Ei and rrii are the energy and mass of the baryon Bi 
respectively. Putting everything together we obtain 



B(B -> A c p) PC = 7.2 x 10" 
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(2.14) 



The PV contribution is expected to be smaller. For example, it is found to be r PV /T PC = 0.59 



in [ r0 1 . Therefore, we conclude that 



B(B° -> A c p) < 1.1 x 10 



-5 



6 



(2.15) 



The strong coupling g B o pS + has been estimated in [I(J using the quark-pair-creation model 
and it is found to lie in the range g B o p ^+ = —(6 ~ 10), recalling that g^NN ~ 14. We 
shall see in Sec. III. A that the measurement of B~ — > A c p7i~ can be used to extract the 
coupling gB+pA b which in turn provides information on g B o pTi +- At any rate, the prediction 
( j2.15| ) is consistent with the current experimental limit 2.1 x 10~ 



21 . Note that all earlier 



predictions based on the QCD sum rule or the pole model [T(J or the diquark model |~T 
are too large compared to experiment (see Table I). In the pole-model calculation in [|10J 
the weak matrix element is largely overestimated. 



^In the bag model the low-lying negative parity baryon states are made of two quarks in the ground 
lSij 2 eigenstate and one quark excited to lP\/ 2 or 1^3/2- Consequently, the evaluation of the PC 
matrix element for ^ — ^ + baryonic transition becomes much involved owing to the presence of 
1Pi/ 2 and I-P3/2 bag states. 
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III. THREE-BODY CHARMFUL BARYONIC DECAYS 



A. B~ — »■ A c pir~ 

The quark diagrams and the corresponding pole diagrams for B~ — > A c p-n~ are shown 
in Fig. 2. There exist two distinct internal W emissions and only one of them is factorizable, 
namely Fig. 2(b). The external W emission diagram Fig. 2(a) is of course factorizable. 
Therefore, unlike the two-body decay B° — > A c p, the three-body mode B~ — > A c p-K~ does 
receive sizable factorizable contributions 

A(B~ - A c p7r-) fact = ^V cb V: d { ai (n-\(du)\0)(A c p\(cb)\B-) 

+ a 2 (n-\(db)\B-)(A c p\(cu)\0)} = A 1 + A 2 , (3.1) 

where naively a± = c\ s + c 2 /N c and a 2 = c 2 +cf^ /N c , to which we will come back later. Un- 
fortunately, in practice we do not know how to evaluate the 3-body hadronic matrix element 

(A c p\(cb)\B~) . Thus we will instead evaluate the corresponding low-lying pole diagrams for 

(*) — 

external lU-emission, namely, the strong process B~ — > A b p, followed by the weak decay 
— > A c 7r~ [see Fig. 2(a)]. Its amplitude is given by 

G 

M = -^V ud V* b g B + P KjnaiU Ac [ft hK \rnl){2p n ■ p Ar + p' 7T (m Ab - m Ac )]75 

+9i" Ac (ml)[2Pn -Pa c -^(m Ab +m Ac )]}v p x ■ — ^ — , (3.2) 

{PAc + Ptt) — 



b 



where we have applied factorization to the weak decay A& — > A c 7r and employed the form 
factors defined by 

(K(PAMcb)\K(PA b )) = U Ac \f^{pl)l, + O^pl + ^^ p^ 

{ f^A b 

-[9^{pl)l» + i^^tf + 9 -^^P^h 5 )u Ab , (3.3) 

where p n = p Ab — p Ac . Note that the \ intermediate state A* h makes no contribution as 
the matrix element (A c |(c6)|A£) vanishes. Likewise, the intermediate states E° and E°* also 
do not contribute to A 1 under the factorization approximation because the weak transition 
(A c |(c6)|E° ( *' ) ) is prohibited as E b and E£ are sextet bottom baryons whereas A c is a anti- 
triplet charmed baryon. 

To evaluate the fcatorizable amplitude A 2 , as shown in Fig. 2(b), we apply the 
parametrization for the B — n matrix element 

(n-(p n )\(db)\B-(p B )) = Ff*tf)(p B + (i^V) - F^(q 2 )) m *~ m * g„, (3.4) 
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FIG. 2. Some pole diagrams for B — > A c pir~ where the symbol • denotes the weak vertex, (a) 
and (b) correspond to factorizable external and internal W-emission contributions, respectively, 
while (c) to nonfactorizable internal P^-emission diagrams. There are two pole diagrams corre- 
sponding to the quark diagram in (c). 



8 



and obtain 



where 



G F 

M = —j=V ud V* b a 2 u Ac [{apn + b) - (c^ + d)j 5 ] v p , 



(3.5) 



a 
b 



2f*"(t)F»*(t) + 2f^(t)F^(t)(m Ac + m p )/m Ac , 



c-A c p 



(m Ac -m p )ft P {t) 



Bit/ 



1 1 

m% — 



-2ft P {t)F^{t){ PAc - p p ) ■ pJm Ac + ft P (t)F Q B *(t)(m 2 B ~ ml)/m Ac , 
c = 2g^{t)F^{t) + 2gt P (t)F^(t)(m Ac - m p )/m Ac , 



d = (m Ac + m p )g^ cP (t) 



Frit) + (*r (*) - F»*{t)) 



m 



B 



m; 



-2gt P (t)Fr(t)(p Ac - p p ) ■ p w /m Ac + <£* {t)F^{t){rn E - m*)/m A 



A c p, 



(3.6) 



and t = q 2 = (p B - p n ) 2 = (p Ac + p p ) 2 . 

The form factors fi and g>j for the heavy-to-heavy and heavy-to-light baryonic transitions 
at zero recoil have been computed using the non-relativistic quark model ||26|| . In principle, 
HQET puts some constraints on these form factors. However, it is clear that HQET is not 
adequate for our purposes: the predictive power of HQET for the baryon form factors at 
order 1/mq is limited only to the antitriplet-to-antitriplet heavy baryonic transition. Hence, 
we will follow |26| to apply the nonrelativistic quark model to evaluate the weak current- 
induced baryon form factors at zero recoil in the rest frame of the heavy parent baryon, 
where the quark model is most trustworthy. This quark model approach has the merit that 
it is applicable to heavy-to-heavy and heavy-to-light baryonic transitions at maximum q 2 and 
that it becomes meaningful to consider l/m q corrections as long as the recoil momentum is 
smaller than the m q scale. It has been shown in [^] that the quark model predictions agree 
with HQET for the antitriplet-to-antitriplet (e.g., A& — > A c , — > H c ) form factors to order 
1/rriQ. For sextet £5 — ► S c and Q b — > Q c transitions, the quark- model results are also in 
accord with the HQET predictions (for details see [p7| 



Numerically we have [27 



fi 



AbA c 



f A b A c/ J2 
J3 



Ai,A c 

9i 

Ai,A c 
92 



ill 



1.02, 
-0.03. 



A b A c 
93 



-0.23, 



(3.7) 



for the A b — A c transition at zero recoil q 2 ^ = {m Ab — m Ac ) 2 , and |26| 



fi cP (<lm 



9i cP {q 2 m 

A c p 
92 



0.80, 
-0.07. 



9s cP (q 2 J = -0.21, 



(3.8) 



for the A c — p transition at q m = (m Ac — m p ) . 

Since the calculation for the q 2 dependence of form factors is beyond the scope of the non- 
relativistic quark model, we will follow the conventional practice to assume a pole dominance 
for the form-factor q 2 behavior: 
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/(* 2 )=w^fgv, 9( rt=.(A)(^^y. (3.9) 



where my (m^) is the pole mass of the vector (axial-vector) meson with the same quantum 
number as the current under consideration. The function 



2 



plays the role of the baryon Isgur-Wise function ((u) for the Aq — > Aq/ transition, namely, 
G = 1 at q 2 = q^. The function ((u) has been calculated in the literature in various different 



models [28-33]. Using the pole masses m v = 6.34 GeV, = 6.73 GeV for the A b — > A c 



transition, it is found that G(q 2 ) is consistent with the earlier soliton model [28] and MIT 
bag model |29| calculation of ((u) for n = 2 PS| . However, a recent calculation of ((u) in 
yields 

CM = (^) (3.1D 

and this favors n — 1. Therefore, whether the g 2 dependence is monopole or dipole for heavy- 
to-heavy transitions is not clear. Hence we shall use both monopole and dipole dependence 
in ensuing calculations. Moreover, one should bear in mind that the q 2 behavior of form 
factors is probably more complicated and it is likely that a simple pole dominance only 
applies to a certain q 2 region, especially for the heavy-to-light transition. For the A c — p 
transition, we will use the pole masses m v = 2.01 GeV and = 2.42 GeV and assume 
dipole q 2 dependence. 

For the form factors F^(q 2 ) we consider the recently proposed Melikhov-Stech (MS) 
model based on the constituent quark picture J53J] . Although the form factor q 2 dependence 
is in general model dependent, it should be stressed that Ff^i^q 2 ) increases with q 2 more 
rapidly than F^ n (q 2 ) as required by heavy quark symmetry. We shall see below that the 
predicted decay rates are insensitive to the choice of form-factor models. 

Thus far we have only discussed factorizable contributions. The nonfactorizable effects 
are conventionally estimated by evaluating the corresponding pole diagrams. The processes 

B - - 7T- + B*° 

^ p + A c 
B~ -> A c + A— 

^P + ti- (3.12) 

are some examples of the pole diagrams shown in Figs. 2(c) and 2(d); they correspond to 
nonfactorizable internal H^-emission. Presumably these nonfactorizable contributions will 
affect the parameter a 2 substantially. 

The total decay rate for the process B~(ps) — > A c (pi) + p{p2) + vr _ (p 3 ) is computed by 
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(2vr) 3 32m| 



1 B 

or 



(2vr) 3 16m| 



|A| 2 dm? 2 dmi 3 , (3.13) 



|A| 2 dS ff dml 3 , (3.14) 



where E 1 ^ is the energy of the outgoing pion, and m?- = (pi + Pj) 2 with p 3 = For a given 
E n , the range of m 23 is fixed by kinematics. Under naive factorization, the parameter ct2 
appearing in Eq. ( |3.1| ) is numerically equal to 0.024, which is very small compared to the 



value of a 2 = 0.40 - 0.55 extracted from 5° -> D°(*h° decays p5f and |a 2 | = 0.26 ± 0.02 



in B — > J/ipK decay As stated before, a 2 may receive sizable contributions from 

nonfactorizable pole diagrams Figs. 2(c) and 2(d). Therefore, we will treat a 2 as a free 
parameter and take a 2 = 0.30 as an illustration. For strong coupling constants a simple 
quark-pair-creation model yields (see Appendix C of |TD| for detail) 



\9B+ P A b \ = 3^3/2 |0 BOpE +|. (3.15) 

Hence, the strong coupling constant \gB+ P A b \ is much larger than \g B o p ^+ |- Putting everything 
together we obtain numerically 

. v ( (lO.Or 2 + 0.04 - 0.8r) x HT 4 for n = 1 

B(B — > A c pn ) — i ; _ / . 

V ; \ (5.5r 2 + 0.04 - 0.6r) x 10" 4 for n = 2 

= f 9.2 x 10~ 4 for n = 1 and 5f B + P A b = 16 . , 

\ 4.9 x 10^ 4 for n = 2 and ^ B+pA() = 16 ' 1 ' J 

where r = g£+ p A b /16 and the first two lines show explicitly the contributions from external 
W-emission, internal jy-emission and their interference, respectively. We find that the ex- 
ternal IV-emission and internal jy-emission contribute destructively (constructively) if the 
A c — p baryonic form factor q 2 dependence is of the dipole (monopole) form. From Eq. ( p. 16 ) 



we find that the strong coupling constant gB+ P A b m the vicinity of order 16 can accommodate 
the observed branching ratio of B~ — > A c p7T~ [see Eq. ( |OJ )]. It follows from Eq. ( p,15| ) that 
\g B o p Y:+ \ ~ 4.3, which is close to the model estimate of 6 ~ 10 given in flOfl . It is likely that 
the quark-pair-creation-model calculation of strong couplings is more reliable for their ratios 
than their absolute values. 

We have checked explicitly that the results are fairly insensitive to the choice of B — 71 
form factors. For example, we have computed the branching ratios using the three different 



form-factor models given in [37] and found that the difference in rates is at most at the level 
of 5%. 

Evidently, the calculated branching ratios are in agreement with experiment ( |1.1| ). There 
are several reasons why the three-body decay rate of B~ — > A c pir~ is larger than that of the 
two-body one 5° — ► A c p. (i) The former decay receives external and internal iy-emission 
contributions, whereas the color-suppressed factorizable IV-exchange contribution to the 
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latter is greatly suppressed, (ii) At the pole-diagram level, the £& propagator in the pole 
amplitude for the latter is of order l/(ml — m 2 ), while the invariant mass of the (A c p) system 
can be large enough in the former decay so that its propagator in the pole diagram is not 
subject to the same 1/m 2 suppression, (iii) The strong coupling constant for A5 — > B~p is 
larger than that for — > ~B^p. 



B. B~ -» A c pp~ 

Naively it is expected that A c pp~ has a larger rate than A c pn~ due to the three polar- 
ization states for the p meson. The calculation for B~ — ► A c pp~ is the same as that for 
B~ — > A c piT~ except that two of the matrix elements are replaced by 

(p-\(du)\0) = f p m p e* p , (3.17) 

and 



(p-\(db)\B-) = e^e^fitfyBptf) - i\ (m B + m fi )e%A?"(q 

Trig ~\- Tflp 

e * 'PB / \ aB P , 2\ n £ ^ 'PS 



2^ 



(p B +P P )»A* P (q 2 )-2m;—^ qil A^(q') - A^(q 2 ) , (3.18 



tub + m p " * r "~ ' '' q 2 



where q = p B - p P = Pa c + p p and 



A B P{q2) = m B + m P A B P{q2) _ m B-m P A B P{q2) _ 



2m p ' 2m p 



Obviously the calculation is much more involved owing to the presence of the four form 
factors V, A , Ax, A 2 compared to the pion case where there are only two form factors F 
and Fx. 

A straightforward but tedious calculation yields 



B(B-^A c pp-) 



(2.6r 2 + 0.02 - 0.3r) x 10~ 3 for n = 1 

(1.5r 2 + 0.02 - 0.2r) x 10~ 3 for n = 2 

2.3 x 10~ 3 for n = 1 and g B + P A b = 16 , . 

1.3 x 10- 3 for n = 2 and g B + P A b = 16 ' [ ' ' 



where we have used the decay constant f p — 216 MeV and the MS model [3H] for the B — p 
form factors. As in the previous case, the contributions from external jy-emission, internal 
iy-emission and their interference are shown explicitly in the first two lines of the above 
equation. Again we have checked explicitly that the predictions are insensitive to the form- 
factor models for the B — p transition. Note that the predicted branching ratio is consistent 
with the current limit on B~ — > A c p7i^7i° [see Eq. (|1.1|) 1. The ratio 
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for n = 1 or 2 is independent of the strong coupling gB+ P A c and hence its prediction should 
be more trustworthy. Experimentally it is important to search for the B decay into A c pp~ 
and have a refined measurement of A c p7r~ in order to understand their underlying decay 
mechanism. 

IV. CONCLUSIONS 

We have studied the two-body and three-body charmful baryonic B decays: B° — > A c p 
and B~ — > A c p7r~(p~). The factorizable W^-exchange contribution to 5° — > A c p is 
negligible. Applying the bag model to evaluate the weak £5 — A c transition, we find 
B(B° — > A c p) < 1.1 x lCr 5 |g£ 0pS +/6| 2 with g B o pTi + being a strong coupling for the decay 

— > and the predicted branching ratio is well below the current experimental limit 
2.1 x 1CT 4 . Contrary to the two-body mode, the three-body decay B~ — > A c pir~ receives 
factorizable external and internal W-emission contributions. The external ^-emission am- 
plitude involves a three-body hadronic matrix element that cannot be evaluated directly. 
Instead we consider the corresponding pole diagram that mimics the external H^-emission at 
the quark level. It is found that the factorizable contributions to B~ — > A c p7i~ can account 
for the observed branching ratio of order 6 x 10~ 4 . The strong coupling \gs+pK b \ is extracted 
to be of order 16, which in turn implies |<7bo P £+| ~ 4.3 under the quark-pair-creation model 
assumption. The decay rate of B~ — > A c pp~ is larger than that of B~ — > A c pix~ by a factor 
of 2.6. We have shown and explained why the 3-body charmful baryonic B decay in general 
has a larger rate than the 2-body one. 

Finally, our present study is ready to generalize to other charmful baryonic B decays, e.g. 
B — > A C A, S c iV, B — > A c A7r(p), T lc Nn(p), ■ ■ ■ , etc. Experimentally it would be interesting 
and important to measure these hadronic decays. 
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APPENDIX 



In this Appendix we evaluate the baryon matrix elements in the MIT bag model |23|. In 
this model the quark spatial wave function is given by 



_ N-i ( ij (xr/R)x 

_ ( iu(r)x \ 
\v(r)a-fx) 



(Al) 



for the quark in the ground (IS1/2) state, where j and j\ are spherical Bessel functions. The 
normalization factor reads 

N-l = 7 7 : TT7^ , (A2) 

[2lo(uj- 1) + mR} 1 / 2 sin x K ' 

where e = (u — mR)/(uj + mi?), x = (lo 2 — m 2 R 2 ) 1 / 2 for a quark of mass m existing within 
a bag of radius R in mode u>. For convenience, we have dropped in Eq. (|A2|) the subscript 
— 1 for x, u and R. The eigenvalue x is determined by the transcendental equation 

x 

taI1X= l_ mjR _( a; 2 + m 2 jR2) l/2- ( A3 ) 

In terms of the large and small components uir) and v(r) of the 1Si/ 2 quark wave function, 
the matrix elements of the two-quark operators V^(x) = q'^^q and A^(x) = q'^^q are given 
by 

(<?'|^o|?) = u'u + v'v, 
(q'\A \q) = —i{u'v — v'u)a ■ r, 
(q'\V\q) = —{u'v + v'u)a x f — i(u'v — v'u)f, 

(q'\A\q) — (u'u — v'v) a + 2v'v fa ■ f. (A4) 

The four-quark operators 0\ = (cb)(du) and 2 = (cu)(db) can be written as 0\(x) = 
6(cb)i(du) 2 and 2 = 6(cu)i(db) 2 , where the subscript i on the r.h.s. of 0\ 2 indicates that 
the quark operator acts only on the 2th quark in the baryon wave function. It follows from 
Eq. (|A4|) that the PC matrix elements have the form 

/ r 2 dr{q' 1 q' 2 \(cb) 1 (du) 2 \q 1 q 2 ) = (-X x + X 2 ) - -pd + ^X 2 )a x ■ <x 2 , 
/ r 2 rfr(g^ 2 |(J6) 1 (c M ) 2 |g 1 g 2 ) = (X, + X 2 ) - l(-X 1 + 3X 2 )a 1 ■ ff 2 , (A5) 
where the bag integrals are defined in Eq. ( |2.11| ) and use has been made of 



dQfifj = -f I dtt, (A6) 
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and those terms odd in r have been dropped since they vanish after spatial integration. Note 
that we have applied the isospin symmetry on the quark wave functions, namely u d {r) = 
u u (r) and v d (r) = v u (r), to derive Eq. (|A5|) . 

We also need the spin-flavor wave functions of the baryons involved such as 

A+ T = -^K cud ~ cdu )XA + (12) + (13)], 

S+T =i&im X . + (12) + (13)], (A7) 
p+T =id««x.+ (12) + (13)], 

where abcXs — (2aWc^ — a^b^c^ — a^b^c^)/y/E, abc\A = (a^b^c^ — a)b^rf)/ y2, and (ij) means 
permutation for the quark in place i with the quark in place j. Applying Eq. ( |A7| ) yields 



(Ac T \bicbibbi d b 2u \^i T) = 0, 



(A+ Ub\ c bubl d h u (a 1 -a 2 )\E+ |) 



1 

76' 



(A c + T \b\MublMK T) = A, (A* 



(A c + T {blhublMia, ■ a 2 )\Et |) 



1 



2y/E' 



In the above equation b\ ql {pi q ) denotes a quark creation (destruction) operator acting on the 
first quark in the baryon wave function. It is easily seen that (A+|Oi|£^) = — (A+|0 2 |£ft~) 
and hence (A+|0 + |E^) = 0, as it should be. The PC matrix element in Eq. (|2.10|) 

<A+|0 PC |£+> = -±(X t + 3X 2 )(4vr) (A9) 

then follows from ( |A5| ) and (|A5|). 

For numerical estimates of the bag integrals Xi and X 2 , we use the bag parameters 

m u = m d = 0, m s = 0.279 GeV, m c = 1.551 GeV, m b = 5.0 GeV, 

x u = 2.043, x s = 2.488, x c = 2.948 , x b = 3.079 , # = 5.0 GeV" 1 . (A10) 



15 



[1] 

[2] 
[3] 
[4] 

[5] 
[6] 
[7 
[8 
[9 
[10 

[11 
[12 

[13 
[14 
[15 
[16 
[17 
[18 
[19 
[20 

[21 
[22 



REFERENCES 

ARGUS Collaboration, H. Albrecht, Phys. Lett. B209, 119 (1988). 

N.G. Deshpande, J. Trampetic, and A. Soni, Mod. Phys. Lett. A3, 749 (1988). 

N. Paver and Riazuddin, Phys. Lett. B201, 279 (1988). 

M. Gronau and J.L. Rosner, Phys. Rev. D37, 688 (1988); G. Eilam, M. Gronau, and 
J.L. Rosner, Phys. Rev. D39, 819 (1989). 

J.G. Kdrner, Z. Phys. C43, 165 (1989). 

A.V. Dobrovol'skaya and A.B. Kaidalov, JETP Lett. 49, 25 (1989). 

V. Chernyak and I. Zhitnitsky, Nucl. Phys. B345, 137 (1990). 

X.G. He, B.H.J. McKellar, and D-d. Wu, Phys. Rev. D41, 2141 (1990). 

G. Lu, X. Xue, and J. Liu, Phys. Lett. B259, 169 (1991). 

M. Jarfi et al., Phys. Rev. D43, 1599 (1991); Phys. Lett. B237, 513 (1990). 

P. Ball and H.G. Dosch, Z. Phys. C51, 445 (1991). 

S.M. Sheikholeslami and M.P Khanna, Phys. Rev. D44, 770 (1991); S.M. Sheik- 
holeslami, G.K. Sindana, and M.P. Khanna, Int. J. Mod. Phys. A7, 1111 (1992). 

G. Kaur and M.P. Khanna, Phys. Rev. D46, 466 (1992). 

CLEO Collaboration, T.E. Coan et al, Phys. Rev. D59, 111101 (1999). 

Belle Collaboration, K. Abe et al, BELLE-CONF-0116. 

I. Dunietz, Phys. Rev. D58, 094010 (1998). 

CLEO Collaboration, S. Anderson et al, Phys. Rev. Lett. 86, 2732 (2001). 
C.K. Chua, W.S. Hou, and S.Y. Tsai, |hep-ph/0107TTO| . 



C.K. Chua, W.S. Hou, and S.Y. Tsai, |hep-ph/010806^ . 



W.S. Hou and A. Soni, Phys. Rev. Lett. 86, 4247 (2001). 

Particle Data Group, D.E. Groom et al, Eur. Phys. J, C15, 1 (2000). 

Y.H. Chen, H.Y. Cheng, B. Tseng, and K.C. Yang, Phys. Rev. D60, 094014 (1999); A. 
Ali, G. Kramer, and CD. Lii, Phys. Rev. D58, 094009 (1998). 



[23] A. Chodos, R.L. Jaffe, K. Johnson, and C.B. Thorn, Phys. Rev. D10, 2599 (1974); T. 
DeGrand, R.L. Jaffe, K. Johnson, and J. Kisis, ibid 12, 2060 (1975). 



16 



[24] H.Y. Cheng and B. Tseng, Phys. Rev. D46, 1042 (1992); 55, 1697(E) (1997). 
[25] H.Y. Cheng and B. Tseng, Phys. Rev. D48, 4188 (1993). 
[26] H.Y. Cheng and B. Tseng, Phys. Rev. D53, 1457 (1996). 
[27] H.Y. Cheng, Phys. Rev. D56, 2799 (1997). 

[28] E. Jenkins, A. Manohar, and M.B. Wise, Nucl. Phys. B396, 38 (1993). 
[29] M. Sadzikowski and K. Zalewski, Z. Phys. C59, 677 (1993). 

[30] X.-H. Guo and P. Kroll, Z. Phys. C59, 567 (1993); X.-H. Guo and T. Muta, Mod. Phys. 
Lett. All, 1523 (1996); Phys. Rev. D54, 4629 (1996). 

[31] Y.B. Dai, C.S. Huang, M.Q. Huang, and C. Liu, Phys. Lett. B387, 379 (1996). 

[32] M.A. Ivanov, V.E. Lyubovitskij, J.G. Kdrner, and P. Kroll, Phys. Rev. D56, 348 (1997). 

[33] M.A. Ivanov, J.G. Korner, V.E. Lyubovitskij, and A.G. Rusetsky, Phys. Rev. D59, 
074016 (1999). 

[34] D. Melikhov and B. Stech, Phys. Rev. D62, 014006 (2001). 
[35] H.Y. Cheng, |hep-ph/0108096| . 



[36] H.Y. Cheng and K.C. Yang, Phys. Rev. D59, 092004 (1999). 

[37] M. Wirbel, B. Stech, and M. Bauer, Z. Phys. C29, 637 (1985); M. Neubert, V. Rieckert, 
B. Stech, and Q.P. Xu, in Heavy Flavours, 1st edition, edited by A.J. Buras and M. 
Lindner (World Scientific, Singapore, 1992), p. 286. 



17 



